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ABSTRACT 

We calculate the evolution of zero-metallicity Population III (Pop III) stars whose mass grows 
from the initial mass of ^ ^Mq by accreting the surrounding gases. Our calculations cover a 
whole evolutionary stages from the pre-main sequence, via various nuclear burning stages, through 
the final core collapse or pair-creation instability phases. We adopt the following stellar mass- 
dependent accretion rates which are derived from cosmological simulations of early structure 
formation based on the low mass dark matter halos at redshifts z ~ 20: (1) the accretion 
rates for the first generation (Pop III.l) stars and (2) the rates for zero-metallicity but the 
second generation (Pop III. 2) stars which are affected by radiation from the Pop III.l stars. 
For comparison, we also study the evolution with the mass-dependent accretion rates which are 
affected by radiatibe feedback. We show that the final mass of Pop III.l stars can be as large 
as ~ IOOOM0, beyond the mass range (140 — 300Mq) for the pair-instability supernovae. Such 
massive stars undergo core-collapse to form intermediate-mass black holes, which may be the 
seeds for merger trees to supermassive black holes. On the other hand. Pop III. 2 stars become 
less massive (< 40 — 6OM0), being in the mass range of ordinary iron core-collapse stars. Such 
stars explode and eject heavy elements to contribute to chemical enrichment of the early universe 
as observed in the abundance patterns of extremely metal-poor stars in the Galactic halo. 

Subject headings: accretion, accretion disks - nuclear reactions, nucleosynthesis, abundances - stars: 
abundances - stars: evolution - stars: formation - supernovae: general 



1. INTRODUCTION 

Just after the Big Bang, a cosmic primordial 
gas consists mostly of H, He and a small amount 
of light elements (Li, Be, B, etc). The first heav- 
ier elements, such as C, O, Ne, Mg, Si and Fe, 
must be synthesized during the evolution of the 



first (metal-free — Population III = Pop III) stars 
early in the history of the universe. It is generally 
thought that massive Pop III stars distribute syn- 
thesi zed heavy elements by supernova explosions 



(e.g., iNomoto et allbood ) 
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If Pop III stars are sufficiantly massive, they 
could also form intermediate mass black holes 
(IMBHs) (- 5 X lO^-^M©). The formation of 
IMBHs would have important implications. Stel- 
lar mass black holes (^ IOM0) are formed as 
the central compact remnants of ordinary mas- 
sive (25 - I4OM0) stars at the end of their evolu- 
tion. Supermassive black holes (SMBHs) (~ 10^ 
- IO^Mq) are now known to exist in the center of 
almost all galaxies (e.g., Kormendv fc Richstond 
1995; Bende^llQOi), but their formation processes 
are largely unknown. Recently, the role of IMBHs 
caught much attention in the context of the for- 
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mation of SMBHs. One of the viable scenarios for 
SMBH formation is a merger tree model in which 
seed black holes of a few xlO^M© formed early 
(e.g., at z ~ 20) are assumed to have merged and 
grown to become SMBHs (e.g., Madau and Reed 



200 ll: IVolonteri et al.M2003l) . In this model, the 



small seed blackholes generally go through IMBH 
stages at some intermediate epochs, e.g., between 
2 ~ 20 and 0. 

For both chemical enrichment and the forma- 
tion of IMBHs in the early universe, the mas- 
sive Pop III stars may play important roles. In 
order to clarify the fate of Pop III stars which 
depends sensitively on the stellar mass, Af , we 
study the evolution of mass accreting Pop III 
stars. In the present paper, we generally call 
stars with mass M > IOOMq very massive stars 
(VMSs). Among VMSs, we refer to stars with 
M > SOO Mpi as core-collapse very massive stars 
(CVMSs) (|Ohkuboet al.ll2006l) . 

The standard cosmological model based on dark 
energy and cold dark matter (CDM), the so-called 
ACDM model, posits that structure forms hierar- 
chically from smaller mass objects (a bottom-up 



scenario 



Davis et al. 1985: Ostriker &: Steinhardt 



200 3t iKirshneiibOoi ). Small density fluctuations 



in the early universe seed nonlinear growth of 
structure via gravitational instability. According 
to this model, smaller, stellar size objects formed 
first and then larger structures such as galaxies 
are assembled by merging of these smaller units. 
The first generation stars are predicted to be 
formed when the age of the un iverse was less than 
a few hunred million yea rs ( Couchman fc Ree^ 
19861 : iTegmark et"al] Il997t lYoshida et al.ll2003h . 
whereas, observationally, t he WMAP data sug- 
gests that it is z > 
Spergel et al]l2007h . 



10 (jSanchez et al. 1 120061: 



In the early universe, dark matter plays an im- 
portant role for the formation of the first gener- 
ation stars. Dark matter (DM) gathers due to 
gravity and forms minihalos with typical mass 
of ~ lO'^MfT^ (H aim an et al.lll996l: iTegmark et"ar 



19971 : iFullerfc Couchmanll2000l ) . After the forma 



tion of minihalos, the baryon gas gravitationally 
collapses to form a star. Since the primodial gas 
contains no heavy elements, molecular hydrogen 
is the only efficient coolant indispensable for star 
formation (Peebles fc Dicke. 1968t Matsuda et al 



19691 : IPallaetal.lll983[ ) 



Theoretical studies on star formation suggest 
that the initial mass function (IMF) of Pop III 
first stars may be different from the present one, 
possibly domina t ed by v ery massive stars (e.g., 
lOmukai & Nishi' 'l998; Nakamur a fc Umemural 



J999; Abel et al. 2002r. lBromm et al. 1999; Omuk ai fc Palla 
|2003|). A number of numerical simulations of 
the formation of the first stars have been car- 
ried o ut. They are, either one-dimension calcula- 
tions ( Haiman et al ]ll996l: iNakamura fc Umemural 
2002t lAhn fc Shapiro 12007 ) o r three dimensional 



simulations (Abel et al 



1999; Bromm fc Larson 



2000L l2002t iBromm et al 



20031: Fuller fc Couchman 



20001: lYoshida et al.l2006l : IO'Shea fc N~manil2006al lbl: 



Gao et alJl2007l) . These studies generally suggest 
that the typical mass of a primodial gas cloud is 
^ 10^ — 1O^M0. The latest detailed calculations 
have shown that a protostellar core of as small 
as M ^ IO^^Mq is formed at the cente r of a 
primordial gas cloud ( Yoshida et al. 2008f) . Af- 
ter such a small proto-stellar core is formed, the 
gas surrounding it accretes on the protostar and 
the stellar mass increases through the pre-main 
sequence stage. 

A crucial parameter is the mass accretion rate, 
which determines the typical mass of Pop III stars 
(e.g. Tan fc McKee 2004) . C osmological sim- 
ulations (e.g. Abel et a l. 2002; Bromm fc Loebl 
2004 lYoshida et aD 120061 : ICao et all l2007l ) have 



shown that the accretion rate is as high as 
M ~ 1O~^M0 yr"-'^ when the core mass is small 
(M < IOMq) and it decreases with increasing 
mass. These studies have suggested that first stars 
might be more massive t han lOOMf:^, in the zero - 
metallicity environment. Omukai fc Palla ( 20031 ) 
calculated protostar evolution with constant ac- 
cretion rates. They treated the accretion rate as 
a free parameter and set it to be of the order of 
M ~ 10"^Mo yr"i and obtained that the final 
mass could exceed 100 Af0. However, they also 
found that, if M > 4 x lO'^M© yr"\ the stellar 
radius rapidly expands to prevent further growth. 
Thus the final fate could not be predicted. If mass 
accretion continues through its lifetime (~ a few 
xlO^ years), i.e., if mass accretion is not impeded 
by feedback from the star itself, the final mass 
may reach several times IO^Mq or more. 

Stars formed in the way d escribed above (which 
are called 'Pop III.l' stars; iMcKee fc Tall l2008h 
radiate a large amount of UV photons and build 
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up HII regions wit h a few kiloparsec diameter 



(jYoshida et al.ll2007n . This environment promotes 
the formation of HD molecules, which can cool the 
gas down to as low as 40 — 50K. In a gas cloud at 
such low temperatures, typical stellar mass formed 
is ~ AOMq, smaller than Pop III.l stars. Stars 
forme d in such environment are called 'Pop III . 2 ' 
stars (jMcKee fc Tanll2008Hjohnson et al.ll2008h . 

The possibility of the formation of VMSs as 
Pop III.l stars has renewed the interest in the 
final fate of such VMSs. Previous studies, which 
did not take into account mass accretion, have 
shown that the final fate of the VMSs is sen- 
sitive to the stellar mass M ( e.s.. iBarkat et al 



Bond et al.|[l984l : lArnettlll996h . If M > SOOM, 



19671: iRakavv fc Shavivl 119681: lOber et all Il983t 



^0, 

the stars undergo core-collapse to form IMBHs, 
which we call CVMS. If M ~ 140 - 300 A/q, the 
stars undergo pair instability supernovac (PISNe) 
and disrupt completely, ejecting a large amount 
of heavy elements. Thus PISNe have been sug- 
gested to be the main source of chemical enrich- 
ment in the early universe. However, recent de- 
tailed comparisons between the ob servations of 



extre mely metal poor (EMP) stars (jCavrel et al 



l2004l) and the nucleosyn t hesis y i elds of PISN mod- 
els (lUmeda fc Nomotol 120021: iHeger fc Wooslev 



20021 ) have shown that the PISN yields are hard to 



reproduce the abu ndance patterns of EMP stars 



(jCavrel et al.ll2004h 



It is interesting to examine theoretically under 
what conditions the Pop III stars end their lives as 
PISNe. The related question of whether CVMSs 
(M - SOO-IO^M©) could actually form is of great 
importance, for instance, to understand the ori- 
gin of IMBHs. For this purpose, in this paper we 
ad opt the realist i c mas s accretion rate obtained 
by lYoshida et al. (|2006 V which follows the evolu- 
tion of dense gas clumps formed at the centers of 
rather low mass (~ lO^^^M©) dark matter halos 
at redshifts z ~ 20. 

In the present paper, we calculate stellar evo- 
lution of Pop III stars growing by accretion and 
investigate the final fate of such stars. So far all 
previous calculations of protostar evolution ended 
at the onset of hydrogen burning, at the start of 
the main sequence, and hence they do not answer 
the question of how such massive Pop HI objects 
evolve after they start hydrogen burning and what 
the final stellar mass can be. Therefore, the ma- 



jor purpose of the current studies is not only to 
follow the pre-main seq uenc e evolution studied b y 
Omukai fc Pallal (|2003h and lYoshida et all (|2006h . 



but also to continue further stellar evolution by 
adopting the mass accre tion rate obtained b y the 
cosmological simulation ( Yoshida et al. 20061 ). We 
calculate whole nuclear burning stages, until the 
point where the star ends its life with core collapse 
and/or explosion. 

Following the introduction in this section, our 
models are described in Section 2 and our results 
are presented in Section 3. In the last two sec- 
tions, 4 and 5, we give discussion, summary and 
concluding remarks. 

2. MODELS IN PRESENT STUDIES 

In Ohkubo et al. (2006), we calculated evo- 
lution, nucleosynthesis, explosion and collapse of 
Pop HI CVMSs starting from the main-sequence 
in the absence of mass accretion. More realis- 
tically Pop HI stars are formed along with cos- 
mological structure formation, and the mass in- 
creases by accretion and the final mass may be in 
the range of VMS (M > lOOM©). In the present 
study, we calculate the evolution of Pop HI stars 
with mass accretion, starting from the pre main- 
sequence phase with small mass comparable to the 
solar mass, until core collapse or pair instabillity 
explosion. In this subsection we summarize our 
models and assumptions at the initial and each 
subsequent stage. 

To calculate presupernova evolution through 
the early hydrodynamical phase, we adopt the 
Heny ey-type h ydrodynamical stellar evolution 
code (iNomoto lll98 2: Nomoto fc Hashimoto 1981 
lUmeda et al.ll999l:IUmeda fc Nomotol2002ll20oi ). 
The mass accr etion is calc ulated wit h the m ethod 
(|l976l ) and iNomoto I (|l982f ). We 



of Neo et al. 



adopt the nuclear r eactio n network developed by 
iHix fc ThielemannI (|l996l ) for calculating nucle- 
osynthesis and energy generation at each stage of 
the evolution. We include 51 isotopes up to Si un- 
til He burning ends, and 240 up to Ge afterwards. 

We start with a l.SM© pre main-sequence star, 
a typical mass size in the low-mass star ranges. We 
investigate how this low mass stellar core grows up 
to a massive one with gas accretion. The starting 
m ass of 1.5M© is larger t han the initial mass set 
bv lOmukai fc Pallal (|2003l ). O.OIM©, by more than 
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two orders of magnitude. However, the time it 
takes for the star to increase its mass from O.OIM0 
to IMq is ~ 10^ yr. This is neghgible compared 
with the overall lifetime of a star. Physically, con- 
vection over the whole star occurs in the contract- 
ing phase. So the starting point little affects the 
later evolution, and as one can see later, we can 
qualitativ ely reproduce th e proto stellar evolution 
shown in lOmukai fc Pallal (|2003h . bmce we con- 
sider Pop III stars the chemical composition cho- 
sen are X(R) = 0.753, X(D) = 2x IQ-^, Xi^Re) = 
2xlO-5,X(4He) = 0.247, and Li) = 2x10-^". 

The most important parameter in this work 
is the mass accretion rate, M. We change this 
value and investigate how the star evolves and 
how large the final mass is. We choose the mass 
accretion rate calculated with three dimensional 
cos mological simul a tions by iYoshida et al. (20(M,) 



and Yoshida et al 



( 20f)7h . These authors calcu- 
lated cosmological structure formation in a ACDM 
universe and evaluated the accretion rate for Pop 
III stars. We approximate M from their results as 
a function of stellar mass M: 



M 



dMy 

dt 



0.0450 X M-2/3 Mgyr-i 

M < 300Mq 
16.3 X M-i-^M0yr^i 

M > 3OOM0 



(1) 



for Pop III.l stars fbv lYoshida et al.ll2006[) . and 



M 



dR'f 



YII 



dt 



0.0250 X M-O-^Moyr^i 

M < IOA/q 
0.10 X M^i-OMoyr-i 

IOA/q < M < AOMq 



M > AOMq 



(2) 

for Pop III.2 stars (bv I Yoshida et al.ll2007[) . These 
accretion rates are shown in Figure [TJ We mark 
the former accretion rate with a subscript 'Y', as 
dMy/dt (Equation [1]). We also adopt models with 
accretion rates smaller than dMy /dt by a factor of 
2, 3, 10, and 20. We call these models 'Y-series'. 
The latter accretion rate (Equation [2]) is labeled 
'Yir, as dMyn/dt (we distinguish this model from 
'Y-series'). For this model, the accretion rate is 
much lower than dMy/dt, because the gas tem- 
perature can become low in the fos sil HII regions 



aroun d the first generation stars (jYoshida et al 
20071 ). In Table [U the models we calculate are 



summarized. 

Protostellar feedback effects may considerably 
affect the mass accretion, hence the protostar's 
growth. Kudritzki (2000) argued that the effect of 
radiation pressure is probably negligible for a Pop 
III star since it has no metal. W e discuss more 
recent studies later in this section. McKee fc Tan 



(|2008[ ) showed that radiative feedback, if signif- 
icant, can stop mass accretion when the stellar 
mass exceeds ~ 100 Af©. If so. Pop III stars around 
this mass range may exist. 

We also adopt accretion rates bv lMcKee fc Tan 
POOSI ). who take interruption by feedback into 
consideration. In their calculation, the accretion 
rate is very high in the proto-star phase, but when 
they reach the main-sequence, the accretion rate 
drops drastically owing to evapolation of the ac- 
cretion disk by ionizing photons. During the main- 
sequence, accretion completely stops when outflow 
and inflow become the same. They parametrized 
a measure of entropy of the accreting gas, and 
subsequent accretion rate varies by more than one 
order of magnitude by changing this parameter. In 
this paper, we adopt typical three accretion rates. 
We ap proximate M from figure 9 in lMcKee fc Tan 
(|2008l ) as follows 



M = 



dMj 



Ml 



dt 



= < 



0.125 X Ar-0-44MQyr~i 

M < SlAfQ 
178 X Af^2.i^^yj.-i 

8IA/0 < Af < 321Af"o 



M > 321Af0, 



(3) 



M = 



and 



dMi 



M2 



dt 



= < 



0.0250 X Af-O-44 ATgyr-i 

Af < 41Af0 
20.0 X Af -2.3 Afoyr^i 

41Af0 < M < 135Af0 



Af > 135Af0, 



(4) 
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M 



M3 



dt 



0.005 X M-0-44^/Qyj.-i 

M < 25Mq 
32 X M-3-2MQyr-i 

25M0 < M < 57M0 



M > 57Mq. 



(5) 

These are marked 'M', as dMui/dt {i — 1,2,3). 
We call these models 'M-series' and we summarize 
them in Table [H These accretion rates are shown 
in Figure [21 

A Pop III star has initially no metal con- 
tent by definition, and thus the line-d riven stel- 
lar w ind is thought to be negligible (jKudritzki 
20001) . There are some mechanisms which lead 
to mass loss - pulsational instability by the e 
mechanism. However, the time scale of amplifi- 
cation of such oscillations is longer than the nu- 
clear burning time scale of the main sequence 
10^ years) for Pop III stars, and so a star 
may evolve wit hout losing a significa.nt fraction of 
its in i tial mass ([Ibrahim et al. 1981 ; Baraffe et al 



20011: iNomoto et all 120031) . On the other hand, 
there are other authors who carried out stellar evo- 
lution with mas s loss even for very lo w-me tallicity 
massive stars. Mevnet et al. ( 2006 ) and Hirschil 
(|2007l ) calculated pre-supernova evolution of very- 



metal poor stars (down to Z = 5 x 10~'^ Zq) with 
rotation and found that such metal-poor stars can 
lose mass because the CNO elements which are 
synthesized in the deep interior are transported 
to the surface due to the rotational mixing. As 
long as mass accretion continues a typical accre- 
tion rate (> lO~'^M0yr~^) is much higher than 
the mass loss rate, and hence in the present study 
mass loss is not considered. 

3. RESULTS 

3.1. Pre Main-Sequence Phase 

Our evolutionary calculation starts from M ~ 
1.5Mq. Figure [3] shows the evolutionary change in 
the stellar radius Ras M increases for each model. 
It shows a few distinct evolutionary phases: pre 
main-sequence phase, main-sequence phase, and 
later phase. Each phase is partitioned with marks. 

Because of the gravitational energy release and 
associated heating, the stellar radius quickly in- 
creases to as large as i? ~ 10^ Rq. At this point. 



M ^ 10~2Moyr~^. After this rapid expansion, 
the star settles into the stage where the radius 
increases only gradually while adjusti ng to the 
mass accretion rate (see also Figure 3 in lNeo et al 



19761 ). At this stage the stellar radius is given ap- 
proximately by the stellar mass and mass accretion 
rate as: 



which is similar to the results b 



(6) 



and Omukai fc Palla (200 



I 



Stabler et al.l 
In the sub- 
sequent evolution, our models exhibit the similar 
trend to those shown in Omukai & Palla (2003; 
Fig.l) and Yoshida et al. (2006; Fig.l3). 

We started our calculations with mass accre- 
tion at M — 1.5Mp,. In mor e realistic calcu l ations 
Omukai k Palial (|2003l ) and lYoshida et all (|2006h 
started their simulation of proto-stellar evolution 
with mass accretion with the very first core mass 
of M ^ lO^'^Af0, although their calculation ends 
at the end of the proto-star phase. The time it 
takes for the star to grow from 10~^Mq to 1.5Mq 
is negligible compared with the lifetime of later 
phases, and this initial stage does not affect how 
massive the star can grow. 

After the proto-star phase, the evolution de- 
pends on the timescales of accretion, Tacc = M/M, 
relative to that of Kelvin-Helmholz (KH) contrac- 
tion, tkh, as shown in Figure |4l 

(1) During the stage with the filled square in Fig- 
ures [3] and m Tacc Tkh, so that the radiative 
luminosity (L) is supplied by the gravitational en- 
ergy release due to mass accretion. The star does 
not contract, i.e., the radius stays almost constant. 

(2) During the stage between the filled square and 
the star mark in Figures [3] and [51 the accretion 
timescale gets longer than the radiative energy 
loss, i.e., Tacc > Tkh- Thus the KH contraction 
releases the gravitational energy at enough rates 
to supply the radiative luminosity. Consequently 
the star contracts as seen as the decrease in the ra- 
dius (Figure[3]). During these stages our results re- 
produ ce a similar trend as found by lYoshida et al 
(|2006h . 



3.2. Main-Sequence Phase 

In the stellar interior the temperature rises dur- 
ing the KH contraction. Because of the absence of 
primordial CNO elements, hydrogen burning takes 
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place through the pp-chain, which does not pro- 
duce high enough nuclear energy generation rate 
in to stop the KH contraction (i.e., Ln < L). 
Eventually, the central temperature (Tc) exceeds 
~ 10^ K and the 3a reaction starts to produce 
"'^^C. When the mass fraction of '^^C reaches 
X(12C) - 10^1" at Tc ~ 10^ K, the CNO cycle 
starts to produce high enough in- The star then 
settles into the main-sequence phase with L ^ L^. 
This is approximately the turning point in the stel- 
lar radius after the KH contraction (see Figure [3]). 
The starting point of the main-sequence is indi- 
cated by the filled star marks in Figure [3l 

For model Y-1, it takes ^ 10^ yr for the star 
to reach the main-sequence, and at this stage the 
stellar mass reaches several x 10 — lOOM©, already 
in massive star range. Model YII is the only case 
in which accretion stops during the KH contrac- 
tion phase. For all other models, accretion con- 
tinues after reaching the main-sequence . Pop HI 
massive stars with M > lOAfQ keep the central 
temperature as high as Tc ^ 1.0 — 1.5 x 10* K 
through the main-sequence. The stellar luminos- 
ity is of the order of ~ 10® — 10^ Lq. The main 
source of nuclear energy generation is the CNO cy- 
cle, and such high temerature is necessary to sup- 
ply high enough nuclear luminosity to keep i ~ in 
with X(CNO) - 10-1° - 10-^ For 'M-series', ac- 
cretion stops during the main-sequence. For YH 
model, accretion stops before the main-sequence. 
After the accretion stops, the star evolves verti- 
cally (i.e., M = constant) in Figure [3l For 'Y- 
series', mass accretion continues throughout the 
evolution. The accretion rate has decreased to 
Tacc ^ ''KH, so that the star changes its structure 
to adjust to the main-sequence structure of the in- 
creased M. Thus the star evolves approximately 
along the main-sequence line in Figure[3]where the 
stellar radius increases with mass {R oc M° ''^). 
The main-sequence phase is the longest in the stel- 
lar life, and thus the final mass is mostly deter- 
mined by how much mass accretes through the 
main sequence. The filled circles in Figure [3] show 
the end of the main sequence (hydrogen exhauston 
at the center) . The star grows to be very-massive. 
The details are sammerized in section 13.31 After 
central hydrogen burning, central helium burning 
follows. The star expands and the stellar radius 
becomes more than 10 times larger at the end of 
helium burning than that on the main sequence. 



Figure [5] shows the evolution of each model on 
the HR diag ram. The se figures correspond to Fig- 
ure 8 in .Omukai fc Pa lla (2003). During the early 
mass accretion phase, the stellar radius gets larger 
and the surface effective temperature decreases to 
as low as ~ 6000K. The H~ bound-free opacity 
is the main source of opacity around this temper- 
ature. After this phase, the star undergoes the 
KH contraction to decrease its radius. The effec- 
tive temperature gets higher to reach ~ 10^ K. 
The main opacity source changes from H" bound- 
free absorption to electron scattering at such high 
temperatures. The star reaches the main-sequence 
and evolves on the HR diagram along the main- 
sequence line with increasing mass. As hydrogen 
is consumed, the star leaves the the main-sequence 
to decrease its effective temperature. At the end 
of hydrogen burning, the stellar mass reaches near 
its final mass. 

Figure [5] shows the evolutionary changes in stel- 
lar lunminosity with increasing mass. For com- 
parison, the main-sequence points of several stars 
with no mass accretion are also shown. After 
reaching the main-sequence, the stars evolve along 
the main-sequence lines of corresponding stars 
with no mass accretion. Figure [7] shows the evo- 
lutionary tracks of the central density (pc) and 
central temperature at later phases through car- 
bon, neon, oxygen, and silicon burning for four 
models. For co mparison, we also show those of 
the 25M0 sta r dUmeda fc Nomotoll2002l ) and the 
IOOOMq star (jOhkubo et al.ll2006l) with no mass 
accretion. Generally, more massive stars have 
higher entropies (lower densities) at the same tem- 
peratures. For the model with mass accretion, 
the central entropy is determined when the stel- 
lar mass is small, and remains low because of 
the small Tacc/TKH even when the stellar mass be- 
comes much larger. After H is ignited, the central 
density-temperature path moves toward higher en- 
tropy, while the central temperature is kept at 
Tc ~ 1 X 10* K. 

3.3. Stellar Lifetime and the Final Mass 

After helium burning (at log Tc (K) ~ 8.5), the 
star with mass accretion evolves almost parallel 
to stars with no mass accretion in the {pc, Tc) dia- 
gram. This is because the evolutionary timescale 
is < 10* yr from central helium exhaustion to core- 
collapse or explosion, and < 1 yr after Tc reaches 



6 



10 K which is much shorter than the accretion 
timescale (race = M/M). Thus, the star evolves 
with neghgible effects of mass accretion, and the 
stellar mass can be regarded as constant, reaching 
its final value M = M{. 

The final mass, Mf, is determined by the stel- 
lar lifetime and mass accretion rate. In Table [2] 
these results are summarized. The final mass for 
'Y-series' is very large. Thi s is primarily due to 
the long lifetime of the star. lOhkubo et al. I (|200d) 
calculated the evolution of CVMS models with 
M = 500Mq and lOOOM© with no accretion. The 
lifetime of CVMS (Af > 300Mq) without accre- 
tion is ^ 2 X 10^ yr, which does not sensitively 
depend on the mass because the M — L relation is 
close to i oc M. On the other hand, our models 
with mass accretion have longer lifetime, > 3 x 10^ 
yr, even if the final mass exceeds 300 Af©. Stars 
with mass accretion evolve through the lower mass 
phase. This makes stellar lifetimes longer. On 
the other hand, for model M-1, the lifetime is as 
short as 2.2 x 10^ yr, much shorter than models 
Y-4 and Y-5, even though the final mass is in the 
same range. For model M-1, accretion stops after 
1 X 10^ yr, just after hydrogen burning starts. At 
this stage the star has already reached the VMS 
range so that the luminosity (and thus hydrogen 
consumption speed) is very high. 

3.4. Final Phases of Evolution and Final 

Fate 

For models Y-1, 2, 3, 4, and M-1, M[ exceeds 
300A//q so that they end their life with iron core- 
collapse. Model Y-5 has Aff = 27?)Mq and ends its 
life as a PISN. For models YII, M-2, and M-3 the 
final mass is too small to become a PISN. These 
stars collapse to form a black hole. Whether a star 
explodes as a PISN or causes iron core collapse is 
decisively important for galactic chemical evolu- 
tion and the origin of elements. Main features of 
the final phases of evolution of these models are 
summarized below and more details will be pub- 
lished elsewhere. 

When central helium is exhausted, the CO core 
contraction accelerates. Figure [8] shows the time 
scale of the increase in the central temperature 
(tt = dt/dlnTc) for models Y-1 (A/f = 915Mq) 
and Yll {M[ = AOMq). From log (K) = 8.6 
to 9.2, Tt decreases drastically from ~ lO'^ yr to 



~ 10"^ yr for model Y-1. At log (K) = 9.2, 
model Y-1 is dynamically collapsing with such 
short Tt, while less massive model YII still un- 
dergoes quasi-static evolution with tt — 10° yr. 
For all cases, the core evolution accelerates, being 
decoupled from the envelope. 

3.4.1. Stars with Mf ~ 140 - 300ArQ 

Model Y-5 ends as a PISN because it enters 
the pair-instability region at log Tc (K) '--^ 9 as 
seen in Figure [71 where the dynamically unstable 
region for the core with the adiabatic exponent 
7 < 4/3 is indicated. The central region collapses, 
which induces thermonuclear runaway of oxygen 
and silicon burning due to the rapid temperature 
increase. The thermal pressure generated by the 



tion and expands the star I 


Rakavv et al. 19671 


Ober et all 1983; Bond et al. 


1984; 


Glatzel et al. 


1985; WooslevI 1986). This appears as the turn- 



ing point of the evolutionary tracks in Figure [71 
The calculation is stopped when the total energy 
of the star becomes positive, which leads to the 
total disruption of the star. 

3.4.2. Stars with Aff > 3OOAf0 

Models Y-1, 2, 3, and 4 enter the pair- 
instability region in Figure [71 However, they do 
not explode as PISNe but undergo core-collapse to 
become black holes. These very massive stars have 
such large gravitational binding energy that the 
nuclear energy release can not exceed this binding 
energy. More details of the evolutionary processes 
of model Y-1 are as follows. 

Figure [9l shows the chemical evolution in the 
interior of model Y-1 (Aff = 915Af0). Helium is 
exhausted in the core around tcoiiapse — t ^ 10"* 
yr. When carbon is ignited at tcoUapse — t ^ 10^^ 
yr, the central region enters the pair instability 
region. 

During collapse (tcoUapsc - t < 10^^ yr), the 
silicon layer and iron core rapidly grow in mass. 
This is because in the intermediate region (Afj. 
100 — 2OOAf0) oxygen and silicon burn explosively 
with the time scale of core collapse. Figure [TOl 
shows the abundance distribution for model Y-1 at 
log Tc (K) = 10.3 and log pc (g cm'^) = 10.0. One 
can see that the mass of the iron core of CVMS is 
large 20 — 25% of the total stellar mass), much 
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larger than ordinary massive stars such as ~ 10 % 
in the 2^Mq star. This r esult agrees wi t h tha t of 



34.3. Stars with Mf < UOMq 

The final mass of model YII (Mf = AOMq) is 
in the range of ordinary massive stars (< 60Mq), 
and so the evolutionary track is similar to that of 
25M0 in the top panel of Figure [7l It ends with 
core-collapse and forms a black hole. 

Model M-2 (Aft = 135Mq) undergoes the core 
oscillation as seen in the {pc,Tc) evolution (Fig- 
ure[71 lower). The evolutionary path in the {pc, Tc) 
plane goes through the marginally stable region 
against the pair-creation instability, so that the 
adiabatic exponent 7 is close enough to induce the 
pulsational nuclear instability due to central oxy- 
gen and silicon burning. In stars with 80Mq < 
Mf < I4OM0, shghtly under the PI SN mass range 



such a core oscillation appe ars (jWooslev et al 



2007t lUmeda fc Nomotol [iooi ) . The nuclear en- 
ergy released during the oscillation is not large 
enough to disrupt the whole star, so that the core 
contracts and central temperature and density rise 
again. The star finally collapses to form a black 
hole after several oscillations (Figure [7]). 

Figure [11] shows the abunance distribution for 
model M-2 (Aft = ISSAf©) at log (K) = 9.95, 
and log Pc (g cm"'^) — 9.3. Contrary to model Y- 
1, the iron core mass is as small as 6AfQ, being 
5% of the stellar mass. One can see from Figure [7] 
that the evolutionary track is similar to that of 
the 25AfQ star at the collapse stage. The silicon 
or oxygen layer does not pass through the pair 
instability region, unlike the VMS star. Neither 
the silicon or oxygen layer can follow the central 
iron core collapse. 

4. DISCUSSION 



In this paper, we have calculated the evolu- 
tion of Pop III stars that undergo mass accre- 
tion. The evolution of Pop III VMS has been 
studied by many authors through later nuclear 
burning stages, collapse, and explosion (e.g. , 
Ober et al.lll983t iBond et al.lll984 lArnettj | l996 : 



Umeda fc Nomotol" 2002t iHeger fc WooslevI 12002 : 
Ohkubo et al.l 2006 ). These studies assumed that 
such stars have their large mass from their start- 
ing point of their life and the same mass has been 



held through the evolution. 

Previous proto-stcUar calculations ended at the 
onset of hydrogen burning on the main-sequence. 
Our results show that with mass accretion the stel- 
lar mass continues to grow after the pre main- 
sequence phase through the main-sequence and 
beyond. Assuming that radiative feedback is not 
significant in the majorit y of cases, w e have con- 
firmed the conclusion of Ohkubo et al.l ([2006) that 
there is a realistic possibility for the existence of 
CVMSs as Pop III stars in the early universe. We 
find that the final mass can be very large, being 
in the range where a star collapses to become an 
I MBH (Aff ^ 300 - 1000 Afg) . 

Ohkubo et all (|2006l) argued that these CVMSs 



contribute significantly to the metal enrichment 
of the early universe, while preventing the over- 
abundance, before ending their life quickly be- 
fore ordinary core-collapse supernovae and hyper - 



Ohkubo et ai] (|2006l) 



novae became dominant 
also found that the abundance of ejected heavy 
metals in this CVMS scenario agrees well with 
the observed chemical abundance data from the 
early universe, in contrast to poor agreement of 
the abundance patterns of PISNe. This is be- 
cause the mass ratio of each element ejected by 
PISNs cannot reproduce the abundance patterns 
of EMP stars, intergalac tic medium (IGM), or in 
tracluster matter (ICM) dUmeda fc Nom oto"2GG: 



Heeer fc WooslevI I2OO2I: Ichieffi fc Limonai 200: 



Dp 

a). 



Furthermore, if the majority of the first genera- 
tion Pop III.l stars are CVMSs, they will form an 
IMBH, not explode as a PISN. Then, the second 
generation (Pop III. 2) stars formed under the in- 
fluence of radiation from Pop III.l stars are less 
massive than ^ IOOATq. This scenario may ex- 
plain why PISN signatures are not observed in the 
abundance patterns. 

The presence of CVMS in the early universe, 
if confirmed, has important implications for pro- 
genitors of IMBHs. There is an interesting sug- 
gestion that some I MBHs may have been fo und 



(jBarth et al1l2005l) . iMatsumoto et all (|200l[ ). by 



using Chandra, reported possible identification of 
a > TOOAfg black hole in M82 as an ultra-luminous 
X-ray source (ULX). More such detections have 
been reported as 
19991 : 



Makishima et al 



Makishimal |2008D 



UL XsdColbert fc Mushotzkv 
2OO0t 



and also as an object in t he Galactic center 



(jHansen fc Milosavlievidl2004D . 
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As to the formation of SMBHs, there are sev- 



eral scenarios (e.g., iReesI I2OO2I 120031 ). SMBHs 
may be formed direct ly from supermassi vc halos 
of dark matter (e.g. , Marchant fc Shapirol Il980t 



In some of these models 
20061 : iTanaka fc HaimanI 



Bromm fc Loebll20Q3tlBegelman et al.ll2006f ) . iMadau 
(|200l[ ) first suggested that mini halos with black 
holes with mass ~ 150Mq at redshift 2 ~ 20 
will merge with each ot her successiv e ly to even- 
tually form SMBHs. Ebi suzaki etal\ (|200l[ ) sug- 
gested a scenario where IMBHs grow to become 
SMBHs by merging and swallowing of many of 
these objects in dense clusters. Several authors 
studied detectability of gravitational waves ra- 
diated by IMBH binaries at high redshifts, by 
observatories suc h as Laser In t erferometer Space 
Ante n na (LISA'idSesana et all l2007t iMicic et al. 



(e.g., IVolonteri fc Rees 
2008h the seed black 
holes are of the order of IOOM0 formed from 
colla pse of Pop HI VMSs, while in some others 



2007 : iTanaka fc HaimanI l2008t IPlowman etal 
2OO9I ). If CVMSs actually existed, they could 



be considered as natural progenitors of IMBHs. 
Our results offer a natural scenario for the forma- 
tion of seed black holes responsible for a merger 
tree model. The existence of CVMS as first stars 
will be supported if in the future neutrinos and 
gravitati onal waves are d etected from collapsing 
CVMSs (jSuwaet al.ll2"009h . 

Our present work uses a result from a three- 
dimensional cosmological simulation by Yoshida 
et al. (2006). Volonteri et al. (2003) have shown 
that the merger tree scenario starting with black 
hole seeds of ~ I5OM0 at z ^ 20 is consistent 
with the present day SMBH population estimated 
from observed quasar luminosity and mass func- 
tions. However, any relevant model must also be 
consistent with the presence of SMBHs massive 
enough (a few x IO^Mq) to power the bright 
high redshift quasars recently discovered by the 
Sloan Digital Sky Survey (see, e.g, a recent re- 
view by Fan 2006). Various models have been 
proposed to explain these high redshift power- 



ful massive SMB Hs (e . g.. IVolonteri fc Reed 



Begelman et al.l 2006t Tanaka fc Haiman 



2006; 



20081 : 



Li et al.l l2007r I2OO8I ). In some of these models. 



the seed black holes have to be formed in hotter 
(virial temperature Tvir > 10'* K) and more mas- 
sive (> IO^Mq) dm halos, and the BHs grow by 
supercritical accret ion, in order to become ma ssive 
enough by z ~ 6. ITanaka fc HaimanI (|2008l ) fur- 



ther extended these earlier studies and examined 
successful scenarios which can be constrained by 
both high redshift and local SMBH populations. 



and R hessceds arc black holes of the order of 10 Mq 
which arc formed directly from DM halos (e.g ., 
Be gelman et al.l 120061: iTanaka fc Haiinaiil I2OO8I ). 
ITanaka fc HaimanI ( 20081 ) conclude that the mod- 
els with seed black holes from both models can be 
successful (to be consistent with both z = 6 and 
local quasar populations), but LISA observations 
can distinguish between them. The environment 
where the seed black holes are formed from more 
massive and hotter DM halos is quite different 
from that of conventional star formation models 
such as those of Yoshida et al. (2006). One possi- 
ble interesting extension of our current work will 
be to adopt supercritical accretion for star forma- 
tion in the environment of more massive hotter 
DM halos. 

Tanaka fc Haiman ( 20081 ) show that in a suc- 



cessful merger tree scenario for z = 6 quasar for- 
mation without supercritical accretion, if the seed 
black holes are IOOM0 formed from Pop HI stars 
and subsequently accrete 60% of the time, 
these holes must be formed rather early at z > 30. 
We note that this constraint may be eased appre- 
ciably if the seed black holes are larger instead, 
such as IMBHs of - lOOOM© from our CVMSs. 
The maximum mass for seed black holes to be 
formed from CVMSs may still become larger if 
the environment of hot, massive DM halos (of 
> 10^ K) is adopted in our future calcula- 



T ■ 

vir 

tions 



5. CONCLUDING REMARKS 

In the present study, we calculate the evolution 
of Pop III stars whose mass grows from the ini- 
tial mass of ^ IM0 by accreting the surrounding 
gases. Our calculations cover a whole evolution- 
ary stages from the pre-main sequence, via various 
nuclear burning stages, through the final core col- 
lapse or pair-creation instability phases. We cal- 
culate models with various mass accretion rates: 
(1) stellar mass dependent accretion rates which 
are derived from cosmological simulations of early 
structure formation (Pop III.l stars), based on the 
low mass dark matter halos at redshifts z ^ 20, 
and (2) mass dependent accretion rates for zero- 
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metallicity but second generation (Pop III. 2) stars 
which are affected by radiation from the first gen- 
eration (Pop III.l) stars. For comparison, we also 
adopt mass dependent accretion rates which are 
affected by radiative feedback. 

The final mass of Pop III.l stars can be very 
large {M ~ IOOOMq), beyond the PISNe mass 
range (14OM0 — 300Mq). Such massive stars form 
IMBHs, which may be the seeds for the forma- 
tion of supcrmassivc black holes. Intriguingly, it 
is theoretically suggested that Pop III. 2 stars are 
less massive (M < 40 — 6OM0), being in the mass 
range of ordinary iron core-coUapsc stars. Such 
stars explode and eject heavy elements to con- 
tribute to chemical enrichment of the early uni- 
verse. The stars in this mass range are favorable 
candidates for the elemental origin of extremely 
metal-poor stars in the Galactic halo. We can ex- 
plain why the signature of PISNe are not seen with 
the senario that Pop III.l stars are very massive, 
i.e., M > 300Mq and Pop III. 2 stars are less mas- 
sive, i.e., M < 40 — 6OM0, although there is some 
uncertainty with radiative feedback. 

The possibility of our scenario that the major- 
ity of Pop III.l stars is CVMSs (M > SOOMq) and 
the majority of Pop III. 2 stars is ordinary massive 
stars [M < 4OM0) is attractive because it explains 
the chemical evolution of the early universe where 
the PISN mass range (140Mq < M < 3OOM0) 
can be avoided. Moreover, the range we found 
for large mass stars offers the attractive possibil- 
ity that these CVMSs are indeed progenitors of 
IMBHs. Although the presence of IMBHs has not 
been firmly established, there are various recent 
observational indications. If many of IMBHs in- 
deed existed in the early universe, even if they 
are rare today, that will provide valuable insight 
to the formation of SMBHs and ultimately galaxy 
formation and evolution. 
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and by the Grant-in-Aid for Scientific Research 
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20674003) and MEXT (17684008, 19047004, 
20040004), Japan. 
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Table 1 

Stellar evolution models with mass accretion rates as a funcstellar mass. The mass 

ACCRETION RATES dMy / dt AND dMyu/dt ARE ILLUSTRATED IN FIGURE [T] THE MASS ACCRETION 

RATES dMui/dt ARE DISPLAYED IN FIGURE [2j 



Models 


Mass accretion rate (MQyi ^) 




Y-1 


dMy/dt 


Equation [T] 


Y-2 


0.5 X dMy/dt 




Y-3 


0.33 X dMy/dt 




Y-4 


0.1 X dMy/dt 




Y-5 


0.05 X dMy/dt 




YII 


dMyu/dt 


Equation [2] 


M-1 


dMui/dt 


Equation [3] 


M-2 


dMM2/dt 


Equation [5] 


M-3 


dMM3/dt 


Equation [5] 



Table 2 

Stellar lifetime, final mass, CO core size, and final fate for each model. 



Models 


lifetime (yr) 


final mass M{ (Mq) 


CO core mass (Mq) 


final fate 


Y-1 


2.3 X 10^ 


915 


370 


Core-CoUpase 


Y-2 


2.4 X 10^ 


710 


330 


Core-CoUpase 


Y-3 


2.5 X 10^ 


610 


275 


Core-CoUpase 


Y-4 


2.9 X 10^ 


385 


170 


Core-Collpase 


Y-5 


3.1 X 10^ 


275 


123 


PISN 


YII 


5.5 X 10" 


40 


14 


Core-Collpase 


M-1 


2.2 X 10'' 


321 


155 


Core-Collpase 


M-2 


3.1 X lO'' 


135 


58 


Core-Collpase 


M-3 


4.5 X lO'^ 


57 


22 


Core-Collpase 



13 




1.5 3 
log M<M^.} 



Fig. 1. — Mass accretion rates for the first generation stars {dMY/dt) and the se cond genera tion stars 
(dMyu/dt) as a function of stellar mass M calculated with cosmological simulations (jYoshida et al.ii2006i) . 
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Fig. 2. — Mas s accre tion rates (dMMi/dt) for the first generation stars as a function of stellar mass M by 
McKee fc Ta3 (|2008l ). The open star marks show the onset of the main-sequence, after which the accretion 
rate drops drastically. The open hexagons show the point where the accretion completely stops. 
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Fig. 3. — Changes in stellar radius, i?, for stars evolving with increasing mass M . The top panel shows 

'Y-series', and the bottom panel shows model YII and 'M-series'. The arrows along the line of model YII 

indicates the direction of evolution. The filled squares on the lines show the beginning of Kelvin-Helmholtz 

(KH) contraction phase (the end of rapid accretion phase) . The filled star marks are the beginning of the 

main sequence. The filled circles show the end of main sequence. 
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Fig. 4. — The evolutionary change in the timescales of accretion, race = M/M, and Kelvin-Helmholz (KH) 
contraction, tkh, for model Y-1 with increasing M . 
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Fig. 5. — The evolution of stars with mass accretion on the HR diagram. Top panel shows 'Y-series', and 
bottom panel shows model YII and 'M-series'. The long dashed li ne shows the location of main sequence stars 
with Z = 0. We adopt the loca tion for M < lOOMg stars from iMarigo et al.l (|200ll ). and for M > lOOM© 



stars from Ohkubo et al. ( 20061 ). 
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Fig. 6. — The evolutionary change in the stellar luminosity with the increasing mass, M, for 'Y-series' (top) 
and model YII and 'M-series' (bottom). The filled star marks indicate the main-sequence of stars without 
mass accretion. 
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Fig. 7. — The evolutionary tracks of the central temperature (Tc) and central density (pc) of models Y-1 
and Y-5 (top panel) and models YII and M-2 (bottom panel). The figures in the brackets are the fina l mass 
Mi. For comparison, the tracks of 25M0 (jUmeda fc Nomotoll200 and IOOOMq (|Ohkubo et al.ll2006f) stars 
are also shown. The dynamically instable regions with the adiabatic exponent 7 < 4/3 are indicated. 
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Fig. 8. — Time scale for central temperature increase, tt, as a lunction of central temperature log Tc after 
helium burning for models Y-1 and YII. 
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9. — Evolution of the chemical structure as a function of time left until the core collapse for model Y 
= 915M0). In each layer the indicated element has the largest mass fraction. 
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